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' New physics in the form of scalar/pseudoscalar operators cannot lower the semileptonic 

the upper bound on the leptonic branching ratio B(B S — > ) sets a strong constraint on 

the maximum value of B(B — > K/j, + [i~) in models with multiple Higgs doublets: with the 

l/"V 

^vq ' current bound, B(B — > K fi + fi ) cannot exceed the standard model prediction by more than 

2.5%. The conclusions hold true even if the new physics couplings are complex. However 
these constraints can be used to restrict new physics couplings only if the theoretical and 
f*) • experimental errors in B{B — > A'/i + ^t~) are reduced to a few per cent. The constraints 

fH ' become relaxed in a general class of models with scalar/pesudoscalar operators. 

■ PACS numbers: 13.20.He, 12.60.-i 
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branching ratio B{B — > K)jl + fi ) below its standard model value. In addition, we show that 



One of the major aims of the large hadron collider (LHC), about to start operating soon, is to 
look for Higgs particles within and beyond the standard model (SM). Even a direct observation of 
a Higgs particle will not suffice to tell us whether it is the SM Higgs or not. An understanding of 

00 . 

possible scalar/pseudoscalar new physics (SPNP) interactions through indirect means is therefore 
extremely crucial. 

The flavor changing neutral interaction b — > sfi + fj,~ serves as an important probe to test higher 



order corrections to the SM as well as to constrain many new physics models. This four-fermion 
interaction is responsible for the purely leptonic decay B s — > for the semileptonic decays 

B — > (K, K*)fi + n~ and also for the radiative leptonic decay B^-^ fj, + fj,~ , y. The semileptonic decays 



have been experimentally observed at BaBar and Belle 
decay has the branching ratio 



,12, 



I, [J] . The pseudoscalar semileptonic 



B(B - A>V~) = (4.2±g;l) x 10" 7 , (1) 
which has been obtained with ~ 350 fb^ 1 of data. These values are consistent with the SM predic- 

jab 

, and the experimental errors are expected to reduce to ~ 2% at the forthcoming 
Super-B factories [u]]. At the moment there is about 20% uncertainty in these SM predictions 
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due to the error in the quark mixing matrix element Vt s and the uncertainties related to strong 
interactions. Improvements in the lattice calculations and the measurement of Vt s are likely to 
bring this error down to a few per cent within the next decade. 

The purely leptonic decay B s — ► is highly suppressed in the SM, the prediction for its 



branching ratio being (3.35 ± 0.32) x 10 -9 Jj|. The uncertainty in the SM prediction is mainly 
due to the uncertainty in the decay constant fp a and Vt s . This decay is yet to be observed in 
experiments. Recently the upper bound on its branching ratio has been improved to [2] 

B{B S -► ^+/0 < 0.58 x 1(T 7 (95% C.L.) , (2) 

which is still more than an order of magnitude away from its SM prediction. The decay B s — > 

will be one of the important rare B decay channels to be studied at the LHC and we expect that 

the sensitivity of about 10~ 9 can be reached in a few years [13] ]. 

In the context of these decays, one needs to focus only on new physics from scalar /pseudoscalar 
interactions, since (i) new physics in the form of vector /axial- vector operators is highly constrained 



by the data on B — ► (K, K*)p + fj, as shown in 



141 ]. and (ii) new physics in the form of tensor 



and magnetic dipole operators does not contribute to B{B S — > A measured value of 

B(B S — ► > 10~ 8 indicates that the new physics must be in the form of scalar /pseudoscalar 

operators. 

We take the effective Lagrangian for the four-fermion transition b — ► to be 

L(b -> = L S m + L SP , (3) 

where 

Lsm = —^VtbV^yjf 1 (s-ff.PLb) fij^/j, + Cio(s7 M P L &) + 2 -^T m b (sia fiu q l/ P R b) fi^fij , 

(4) 

L S P = ^V tb vARs{sP R b)^^ + R P (sP R b)ji lh A . (5) 
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Here Pl,r = (1 =F 7s)/2 and q is the sum of the fi + and pT momenta. Rs and Rp are the scalar 
and pseudoscalar new physics couplings respectively, which in general can be complex. We use 
the notation Rs = Rse t&s ,Rp = Rpe lSp . Here the phases are restricted to be < (8s, <5p) < vr, 
whereas Rs and Rp can take positive as well as negative values. Within SM, the Wilson coefficients 
in Eq. (JH) have the following values (|: 

Cf = -0.310 , Cf = +4.138 + Y(q 2 ) , C w = -4.221, (6) 



where the function Y{q ) is given in [15|]. These coefficients have an uncertainty of about 5%, 
which arises mainly due to their scale dependence. 

In Eq. (J5|), we have taken only Pr in the quark bilinear, while the most general Lagrangian must 
have a linear combination of Pl and Pr. Here we start by considering the simpler case because 
SPNP operators mostly arise due to multiple Higgs doublets. In such models, the coefficient of Pr 
in the Lagrangian is much larger than that of Pl [(J. In two Higgs doublet model, for instance, 
the coefficient of Pl is smaller by a factor of m s /mb [if] ]. We shall examine the consequences of 
considering the most general quark bilinear in the latter part of this Letter. 

In the following, we consider the interrelations between the contributions of L$p to the branching 
ratios of the decays B s — > /x + /U~ and B — ► Kfi + fi~ . The effect of SPNP couplings on additional 
observables related to these decays, viz. forward-backward asymmetry in the semileptonic decay 
and the polarization asymmetry in the leptonic decay, has been studied in [jj]. The contribution 



of Lsp to B — > K*fi + /i is so small [g] that no worthwhile correlation can be established between 
it and other decays. Also, Lsp does not contribute to the radiative leptonic decay B s — > /x + ^~7 

0,Q. 

We first consider the contribution of L$p to the decay rate of B s — > The branching ratio 

is given by 

Bsp(B s -> /,+/*-) = 647r ? - \Vt b V t *\ 2 fl(R 2 s + R 2 P ) • (7) 

Taking f Bs = (0.259 ± 0.027) GeV Q, we get 

B SP (B S -> (i + fT) = (1.43 ± 0.30) x 10~ 7 (R 2 S + Rp) . (8) 

Note that the present experimental upper limit on B(B S — ► fi + [t~) is an order of magnitude 
larger than the SM prediction. In the following, we will assume that the SPNP will provide an 
order of magnitude increase of B(B S — > In such a situation, the SM amplitude can be 

neglected in the calculation of the branching ratio. Equating the expression in Eq. ([HJ to the 
present 95% C.L. upper limit in Eq. ([2|), we get the inequality 

{R 2 S + R 2 P ) < 0.70 , (9) 

where we have taken the 2a lower bound for the coefficient in Eq. (JSj) . Thus, the allowed region 
in the Rs~Rp parameter space is the interior of a "leptonic" circle of radius rg ~ 0.84 centered at 
the origin, as indicated in both the panels of Fig. [TJ As the upper bound on B(B S — > /i + ^~) goes 
down, the radius of the circle will shrink. 
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FIG. 1: The allowed ranges of Rs and Rp, when the new physics couplings are real. In both figures, 
the dark grey circles centered at origin represent the regions allowed by the current 2a upper bound on 
B{B S — > fi + The light grey annulus in each figure represents the parameter space allowed by B(B — > 
K/i + n~) at 2cr. The width of the annulus corresponds to the sum of the theoretical and experimental 
errors, both of which are taken to be 2%. In the left panel, we take B(B — ► K^i + ^i~) = (5.64 ± 0.11) • 10 -7 . 
The overlap between the allowed regions is represented by the black crescent. In the right panel we take 
B(B — > KjjL + n~) = (6.04 ± 0.12) • 10 -7 , where the allowed parameter spaces do not overlap. 



We now turn to the semileptonic decay B — > K[i + n~ . The measured branching ratio is consis- 
tent with the SM prediction, though there is a 25% error in the measurement and about 20% error 
in the theoretical prediction due to uncertainties in Vt s , form factors and Wilson coefficients (which 
in turn depend on Vt s ). With the addition of the SPNP contribution, the theoretical prediction 
for the net branching ratio becomes [f| 

B(B -> Kn+yT) = [5.25 + 0.18(12! + R P )-0.13R P cos 5 P ] (1 ± 0.20) x 1(T 7 , (10) 

In Eq. pop, the first term is purely due to the SM, the second term is purely due to SPNP and 
the third term is due to the interference of the two. The theoretical errors arise from one tensor 
and two vector form factors in the SM, and a scalar form factor in SPNP (which is related to 
one of the SM vector form factors). We have made the simplifying assumption that the fractional 
uncertainties in all the form factors are the same. 
Eq. pop can be rewritten as 

B{B Kfj+fi') = (1 + e)B su , (11) 

where e is the fractional change in the branching ratio due to SPNP. The maximum negative value 
that e can take is —0.005, thus implying that the SPNP new physics cannot lower the branching 
ratio B(B — > Kfi + fi~) by more than 0.5% below its standard model value. Indeed, if the theoretical 
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and experimental errors in this quantity were improved to 5%, with the central values unchanged, 
the discrepancy cannot be accounted for by SPNP at 2a. 

Let us first consider the case where the new couplings Rs and Rp are real, which is typical for 
the class of models where the only charge-parity violation comes from the CKM matrix elements. 
Using Eqs. dJJ and (fTUl) . we get 

R% + (R P - 0.36) 2 = -2*- = - 29.04 , (12) 

s y 1 (0.18 ±0.036) x 10~ 7 v ; 

where B cxp is the measured value of B(B S — > Kfi + fi~). The region in the Rs~Rp plane allowed 
by the measurement of B{B S — ► Kfj, + ^~) is then an "semileptonic" annulus centered at (0,0.36), 
as shown in both the panels of Fig. [TJ The inner and outer boundaries of this region correspond 
to the lower and upper bounds of the right hand side of Eq. (|12|) . The right hand side turns out to 
be negative if B cxp is below the SM prediction by more than 0.5%. Then the radius of the circle 
becomes imaginary, which implies that the discrepancy of the measurement with the SM cannot 
be explained by SPNP. 

To illustrate the tension between the quantities B(B S — > ^ + /i~) and B(B — > Kfi + fi~), we 
consider the scenario where the errors in both Bsm an d -B C xp have been reduced to 2%, while 
keeping the upper limit on B(B S — > fj> + n~) at its current value. The allowed Rg-Rp parameter 
space is shown in Fig. [TJ If the lower limit on B exp is small enough, the semileptonic annulus will 
overlap with leptonic circle, as shown in the left panel. However, if the lower limit on -B cxp is larger 
than a critical value (determined by the bound on the leptonic branching ratio), then there is no 
region of overlap as shown in the right panel. In such a situation, the difference between B exp and 
Bsm cannot be accounted for by SPNP because of the constraint coming from the leptonic mode. 

We represent the radius of the leptonic circle by 77 and the inner (outer) radius of the semilep- 
tonic annulus by rj n (r out ). There is tension between the two measurements if 

r in -r e > 0.36 , (13) 

in which case the regions allowed by the two branching ratios do not overlap. Given the current 
value of ri = 0.84, we require < r^ n < 1.2 for an overlap. This implies that the 2a lower limit 
on B exp should be between 4.93 x 10~ 7 and 5.67 x 10~ 7 . (We have added the theoretical and 
experimental errors in quadrature.) If the upper bound on B(B S — > fi + fj,~) is improved by a factor 
of 5, the 2a range for the lower limit on B cxp would be (4.93 — 5.57) x 10 -7 . For the tension to be 
manifest in future experiments, the reduction of errors in B cxp and Bsm is the most crucial. 



When Rs and Rp are complex, the constraint Eq. (I12j) becomes 

R 2 S + {Rp - 0.36 cos 5p) 2 = ^^gL— — - 29.17 + (0.36 cos 5 P f . (14) 

For nonzero 6p, the center of the semileptonic annulus shifts along the Rp axis, while the radius 
of the annuli are almost unchanged. If the allowed regions do not overlap for 5p = (as illustrated 
in the right panel of Fig. [I]), then they will not overlap for any value of dp. Hence the tension 
between B(B S — ► and B{B — > Kfi + fi~) persists, and gives rise to the same constraints on 

the semileptonic branching ratio even if the SPNP couplings are complex. 

In writing the effective SPNP Lagrangian in Eq. ((^J), we considered only the quark bilinear 
sPpb. Lorentz Invariance of the Lagrangian also allows the bilinear sPzb in general. We can take 
this generalization into account by replacing sPpb by s(aPi + Pp)b, where a is the strength of the 
sPjJ) bilinear relative to that of sPpb. With this modification, B(B — > K[i + fi~) is driven by the 
sum of the two quark bilinears with different chiralities, whereas B(B S — > /i + /i _ ) depends on their 
difference {^(J. The expressions for the branching ratios of the two processes considered here are: 

B(B S -» n+iT) = (1 - a) 2 (R 2 s + Rp) (1.43 ± 0.30) x 10~ 7 , (15) 
B(B — > Kfi + fi~) = [5.25 + 0.18(1 + a) 2 {R 2 S + Rp) - 0.13 (1 + o)22 P ] (1 ± 0.20) x 10~ 7 .(16) 

Here we have taken Rs,Rp and a to be real for simplicity. For a = 0, Eqs. (|15p and (|16|) reduce to 
Eqs. ([8]) and (|10p respectively. For the special case a = 1, the new physics has no contribution to 
B s — ► /-f + /U~ because the quark bilinear is pure scalar and the corresponding pseudoscalar meson 
to vacuum transition matrix element is zero. In such cases, B{B S — > /i + ^~) is entirely due to the 
SM, and provides no constraints on B(B — > Kfi + fj>~). 

In Fig. [21 we show e max , the maximum fractional deviation of £>(i? — > Kfi + fi~) from its SM 
prediction as defined in Eq. (jlip . as a function of the 2a upper bound on B(B S — > fi + fJ,~). The 
minimum allowed value of e is almost independent of the value of a and the leptonic upper bound, 
and is approximately —0.005. For the class of models with multiple Higgs doublets, a = 0, and 
the maximum value of e is restricted to +0.025, as seen in earlier discussions. With the additional 
freedom generated by the extra parameter a, this severe constraint is relaxed. For example, for 
the models with a ~ 1.5, the value of e may be as large as +0.7, as can be seen in the figure. In 
general for positive a values, e max increases with a for a < 1.0, and decreases thereafter. When 
a < 0, Eq. (|15p indicates that the constraints on Rs and Rp should become more restrictive. As 
a result, e is constrained to be even smaller. From the figure, e max for negative a are seen to be 
very close to zero, and the corresponding e max curves are almost overlapping. This implies that for 
negative a, any significant deviation of B(B — > Kfi + fi~) from SM is impossible with SPNP. 
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FIG. 2: e max as a function of the 2a upper bound on B(B S — > for different values of a. We have 

taken the theoretical error on B(B S — > fi + fj>~) to be 20%; decreasing it would further constrain e max . 



For the measurements of B(B S — > and B(B — > Kfi + [i~) to be compatible with SPNP, 

the lower bound on B(B — > K/j, + fi~) should be less than (1 + e max )i?sM- Thus, the upper bound 
on B(B S — > and the lower bound on B(B — > Kfi + fi~) allow us to constrain the value of a 

in a class of models that involve new physics scalar /pseudoscalar couplings. 

In this letter, we have parameterized scalar /pseudoscalar new physics in terms of the effective 
operators given in Eq. ([5]). In general, the introduction of new scalar /pseudoscalar fields into a 
model leads to not only new effective operators but also modification of the coefficients of the SM 
operators, e.g. the Wilson coefficients CV, Cg and C\q shown in Eq. However, it has been 
shown that these modifications due to new scalar /pseudoscalar fields are very small [l^. We 
have computed these changes in the two Higgs doublet model and found them to be at most 1%. 
Thus, our assumption of retaining the SM values for the Wilson coefficients, even in the presence 
of new scalar /pseudoscalar fields, is valid. 

In summary, we have shown that in a class of models with new scalar /pseudoscalar operators, 
which includes models with multiple Higgs doublets, the SPNP couplings are strongly constrained 
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by the upper bound on B(B S — > fi + fj,~), and in turn restrict the allowed values of B(B — > Kfj, + fi~) 
to within a narrow range around its SM prediction. Future precise measurements of these two 
branching ratios have the potential not only to give an evidence for new physics, but also to reveal 
the nature of its Lorentz structure. However in order to achieve this, the theoretical as well as 
experimental errors on B{B — > Kfi + fi~) need to be reduced to a few per cent. 
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